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Abstract 

Non-extensive statistical physics has allowed to generalize mathemat- 
ical functions such as exponential and logarithms. The same framework 
is used to generalize sum and product so that the operations allow a 
more fluid way to work with mathematical expressions emerging from 
non-additive formulation of statistical physics. In this work we employ 
the generalization of the exponential, logarithm and product to obtain a 
formula for the survival fraction corresponding to the application of several 
radiation doses on a living tissue. Also we provide experimental recom- 
mendations to determine the universal characteristics of living tissues in 
interaction with radiation. These results have a potential application in 
radiobiology and radiation oncology. 

Entropy, Statistical mechanics. Radiation, Non-extensivity 

1 Introduction 

Biological systems are the paradigm of open non equilibrium systems charac- 
terized by long range correlations and by constantly occurring non markovian 
processes. This makes its description in terms of equilibrium statistical physics 
a very hard task. Thus a formulation of processes in this kind of systems must 
be made in a non orthodox statistical framework in which the concept of entropy 
is expressed in a more general way. 

Tsallis's definition of entropy has been successfully used, in the last years, 
to extend the classical statistical mechanics to several problems involving very 
complex correlations and properties. Its expression for the entropy (in units of 
the Boltzmann's constant) [1] is 

s,=-^(^i- p''{x)dxy (1) 



1 



where p{x)dx is the probabiHty of the magnitude x to adopt values between 
X and X -\- dx and q is the "index of non-extensivity". This, together with the 
constraints that the probabiHty density satisfies the normaUzation condition and 
the finiteness of the "g-mean value" p''{x)xdx = {x) ^ < oo, has been applied 
in different problems with excellent results even when statistical physics based 
on Boltzmann-Gibbs (BG) entropy had not succeed. On the other hand, Tsallis 
entropy reduces to the BG one when q = 1, so that in this sense the former 
entropy can be interpreted as a generaHzation of the last one. 

Our objective in this work is to employ this new formulation of statistical 
physics and the consequent generalization of several mathematical functions 
and operations, provided in the framework of this new formulation, to obtain a 
simple universal formula for the survival fraction of tissue cells under radiation. 

We apply the maximum entropy principle using the Tsallis entropy to obtain 
an expression for the survival fraction in terms of the radiation dose. Then using 
a generalization of the definition of logarithm and exponential function and a 
generalized expression for the product the resultant survival fraction of several 
radiation sessions is obtained. 



2 Materials and methods 

To apply the maximum entropy principle in the Tsallis version to the problem 
of finding the survival factor of a living tissue [2] that receives a radiation, we 
postulate the existence of some amount of absorbed radiation A < oo{or its 
equivalent "Hmit effect", O = agA) after which no cell survives. The application 
of the maximum entropy principle performs like the usual one but with a few 
modifications. 

The Tsallis entropy becomes: 

S,^-^(^l- j\'^{E)dEj , (2) 

where p{E) is the cell death probability, dependent on the "radiation effect" 
E, which is a function of the applied dose D. The way E depends on D is a 
controversial matter, being generally accepted a linear relation E = aD and 
also the "linear quadratic model" E = aD + . 

The normalization condition is in this case £ piE)dE = 1 and the g-mean 

value becomes p'^{E)EdE — (E)^ < oo. With this definition, all properties 
of the tissue and its characteristics of the interaction with radiation become 
included in (E)^ and therefore in il. This is the only parameter (besides q) 
entering in our description. It is clear that the determination of {E)^ for the 
different tissues under different conditions of radiation would give the necessary 
information for the characterization of the survival factor. 

To calculate the maximum of ^ under the above conditions the well known 
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method of Lagrange multipliers [3] is applied, obtaining 



and 





(3) 



(4) 



Then the survival factor is 

Fs{E)^ j\{x)dx^{l-^''' (5) 

with q <liov E <n and zero otherwise. It is not hard to see that when g — > 1 
then ^ cx) and {E) ^ {E). 
Equation ([5]) can be written 

^ ^ [O VD ^ A ^ ' 

where we introduced E = uqD and 7 = 

Equation ([6]) represents the survival fraction in terms of the measurable 
quantities D (radiation dose) and A. All the information about the kind of 
radiation, radiation rate, etc is contained in the phenomenological term A, 
whereas tissues are characterized by 7. This makes ^ a very general expres- 
sion with universal characteristics since the phase transition described by ([6]) is 
homomorphic with the phase transition of ferromagnets near the Curie point. 
The exponent 7 in this case, as in ferromagnetic phase transitions, determines 
the universality class. Then 7 in our case deals only with the kind of tissue that 
interacts with radiation [4j. 

The value of A characterizes a critical point for cell survival probability 
such that for < A the probabilities of cell survival and death coexist but 
when D > A a phase transition takes place and no cell survives. Note that 
we have postulated a linear relation between E and D. With this, the curious 
property that the tissue effect E is additive while the survival fraction is not 
multiplicative emerges. 

The linear model for the tissue effect j5j, however, gives that if the dose is 
additive the corresponding survival fraction is multiplicative if their mean values 
are the same. Though this property belongs only to the linear model and not to 
more general descriptions Hke the LQ model [5j and others, we think it is worth 
to find a link between the additivity property of the dose and the probabiHstic 
properties of the cell survival fraction. 

To motivate this let us first observe that the Tsallis definition of entropy, 
though gives a more general expression than the BG one, lacks one of the im- 
portant properties of this entropy: extensivity. 



3 



Indeed, let us consider a system composed by two independent subsystems A 
and B. Extensivity of the entropy means that the entropy of the whole system 
is the sum of the entropies of the subsystems: S{A U B) — S{A) + S{B). But 
in the Tsallis case Sq{A U B) = Sq^A) + Sq^B) + {q - l)Sg{A)Sq{B). 

As Tsallis entropy can be defined through the introduction of the g-logarithm: 



Sq — 



p^{x)lq {p{x)) dx 



(7) 



' T^^^ ^ , it is possible to generalize the operation of multiplication 
introducing the g-product ®q such that x®qy = expq [lq{x) + lq{y)] [6J where 



where lq{x) 



expq{x) = [! + (!- q)x\ i-" 



(8) 



Taking into account that for the survival fraction we obtain expression ^ 
its analogy with ^ is evident. Then, let us define the exp^{x) function 



exp-y^x) — 1 + 

and the l^{x) its inverse function: 

l^{expj{x)) — X. 
Let us introduce the 7-product of two numbers x and y as: 



x<?)~fy = expj ihi^) + h{y)] = 



XT + — 1 



(9) 



(10) 



(11) 



Note that definitions ^ and (flOl) are not essentially different from Cq {x) and 
lq{x). We are just introducing these definitions to simplify the calculations. 



3 Results 

Let us now define the "effect potential" as x = —lj{Fs). We demand this 
potential to satisfy the additive property. Then the survival fraction expressed 
as 



Fs = exp-^{~x) = exp^( — 



(12) 



where the effect potential depends on (E)^, becomes 7-multipHcative. This 
impHes that the statistical independence of the survival fractions is only possible 
when 7 ^ 00 {q ^ 1). 

Then the survival fraction for the sum of two potentials i.e., the survival 
fraction of a dose B and a dose A is: 



FsiXB) = 



1 - 



XA + XB 

7 



(13) 
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and for n doses 



Fsinx) = 


n 


7 




7 

= expj 


n 




n 

(8) 














.1=1. 
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where [0-Li]^ denotes the iterated appHcation of the 7-product. 

4 Discussion 

The obtained result permits to calculate directly the result of the application 
of different radiation doses if 7 and A are known. Then the importance of 
experimental determination of 7 and A becomes evident. 

To determine A means to measure the minimal dose for which, in a given 
tissue, no cell survives. This can be made by a progressive application of a given 
kind of radiation. Once A is known in those radiation conditions, a log-log plot 
of Fs versus 1 — D/A must be fitted with the straight line whose slope is 7. 
Once 7 is found the tissue is already characterized. Any other kind of radiations 
or the same radiation applied in different conditions to the same tissue can be 
fitted to the same straight line. 

This way, the introduction of the rules for generalization of the product as in 
[6] permits an easy interpretation of the properties of the survival fraction and a 
formula to calculate it for several radiation doses: the survival fraction of several 
radiation doses is the 7-product of the separate survival fractions corresponding 
to each dose. 
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